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ABSTRACT 

The homeobox transcription factor CDX2 plays a 
crucial role in intestinal cell fate specification, both 
during normal development and in tumorigenic 
processes involving intestinal reprogramming. The 
CDX2 regulatory network is intricate, but it has not 
yet been fully uncovered. Through genome-wide 
screening of a 3D culture system, the RNA-binding 
protein MEX3A was identified as putatively involved 
in CDX2 regulation; therefore, its biological rele- 
vance was addressed by setting up cell-based 
assays together with expression studies in murine 
intestine. We demonstrate here that MEX3A has a 
repressive function by controlling CDX2 levels in 
gastric and colorectal cellular models. This is 
dependent on the interaction with a specific 
binding determinant present in CDX2 mRNA 
3 untranslated region. We have further determined 
that MEX3A impairs intestinal differentiation and 
cellular polarization, affects cell cycle progression 
and promotes increased expression of intestinal 
stem cell markers, namely LGR5, BMI1 and MSI1. 
Finally, we show that MEX3A is expressed in 
mouse intestine, supporting an in vivo context for 



interaction with CDX2 and modulation of stem cell 
properties. Therefore, we describe a novel CDX2 
post-transcriptional regulatory mechanism, 
through the RNA-binding protein MEX3A, with a 
major impact in intestinal differentiation, polarity 
and sternness, likely contributing to intestinal 
homeostasis and carcinogenesis. 

INTRODUCTION 

The homeodomain transcription factor CDX2 is a critical 
determinant of intestinal homeostasis, both during devel- 
opment and throughout adult life. CDX2 is involved in 
the antero-posterior patterning of the mammalian embryo 
and is the key molecular mediator of intestinal differenti- 
ation (1^1). Furthermore, multiple evidences substantiate 
CDX2 crucial role in carcinogenesis of the digestive tract. 
It was shown to inhibit cell growth and migration in vitro 
as well as dissemination of colon tumour cells in vivo (5). 
CDX2 heterogeneous loss has also been observed in colo- 
rectal carcinomas (CRCs), particularly in invasive cells 
at the tumour edge (6). Moreover, CDX2 reduction 
increases the progression of chemically induced CRCs 
(7). Conversely, under certain pathological conditions, 
CDX2 becomes abnormally expressed in other organs of 
the digestive tract besides intestine, namely the esophagus 
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(8) and stomach (9,10), driving a precancerous lesion 
known as intestinal metaplasia, a process confirmed in 
transgenic mouse models (11,12). 

Owing to the essential function in intestinal develop- 
ment, differentiation and carcinogenesis, CDX2 regula- 
tion has been extensively studied. We have previously 
identified different mechanisms involved in the tran- 
scriptional regulation of this gene such as the Bone 
morphogenetic protein (BMP) pathway (13), SOX2 (14) 
and a CDX2 autoregulatory loop (15). Several transcrip- 
tion factors including HNF4a, GATA6, TCF4 and 
P-catenin were shown to interact with Cdx2 promoter 
fragments (16). However, mutations at the CDX2 locus 
are a rare event in CRC (17), and its expression does 
not depend on methylation of the proximal promoter 
(18). On the other hand, CDX2 protein phosphorylation 
has also been shown to modify its activity in intestinal 
cells (19,20). These and other studies support the notion 
that CDX2 regulation is intricate and strictly controlled. 

During the past two decades, post-transcriptional regu- 
lation emerged as a fundamental mechanism guiding gene 
expression in higher eukaryotic cells, being at the core of 
normal cellular processes but also cancer initiation and 
development. It is now increasingly clear that RNA matur- 
ation, localization, translation and stability provide 
multiple layers of spatio-temporal control determining a 
transcript's fate (21,22). These coupled events are generally 
dependent on the cooperation between m-regulatory 
elements and fram-acting factors, such as non-coding 
RNAs and RNA-binding proteins (RBPs). RBPs have 
been implicated in virtually every aspect of RNA metabol- 
ism (22,23), particularly, their repressive role is critical to 
establish precise translational patterns that define develop- 
mental and differentiation switches in many organisms. 

In Caenorhabditis elegans, MEX-3 is a translational 
repressor that regulates blastomere identity during early 
embryogenesis (24) and germline totipotency in the adult 
worm (25). MEX-3 has two K homology domains, which 
are conserved single-stranded RNA-binding motifs (26). 
Mutations disrupting mex-3 locus are lethal, resulting in 
embryos that inappropriately generate body-wall muscle 
from the anterior blastomere; hence, the name mex for 
'muscle excess'. This is specified, in part, by the repressive 
function that MEX-3 exerts over the transcription factor 
pal-1, the C. elegans orthologue of caudal in Drosophila 
and CDX in mammals (24,27). In humans, mex-3 was 
identified and characterized as four homologous genes, 
MEX3A-D (28,29), whose biological relevance is starting 
to be explored. Recently, the functional role of MEX3C 
as a RNA-binding ubiquitin E3 ligase was established, 
mediating the post-transcriptional decay of HLA-A 
allotypes (30). It was also shown to be necessary for 
normal postnatal growth in mutant mice by enhancing 
the local expression of insulin-like growth factor 1 in 
bone (31) and appears to be involved in metabolic regu- 
lation of energy balance (32), through yet unknown 
effectors. A variant form of MEX3D called TINO was 
shown to negatively regulate the antiapoptotic protein 
BCL-2 in HeLa cells (33). Finally, knockdown of 
MEX3A by siRNAs was shown to suppress cell prolifer- 
ation and migration in human gastric cancer cells, but the 



molecular mechanisms behind these findings were not ad- 
dressed (34). Pursuing the aim of uncovering new CDX2 
regulatory mechanisms, we explored a putative transla- 
tional repression by MEX3A, which was inversely 
correlated with CDX2 in a 3D experimental model. By 
studying MEX3A expression in vivo and using a cell 
line-based approach to modulate its levels in vitro, our 
study describes a novel post-transcriptional process by 
which CDX2 expression is impaired in the gastrointestinal 
setting and intestinal-like homeostasis compromised, 
through alterations in differentiation, polarity and 
sternness features. Another layer of control is thus added 
to the complex CDX2 regulatory network, involving 
MEX3A as a key regulator of intestinal homeostasis, 
which might have significant implications to gastrointes- 
tinal carcinogenesis. 

MATERIALS AND METHODS 

Cell culture and treatments 

Human gastric carcinoma cell line AGS (ATCC, 
American Type Culture Collection) and CRC cell line 
Caco-2 (ATCC) were cultured under standard conditions 
in RPMI-1640 medium and Dulbecco's modified Eagle's 
medium, respectively, containing 10% fetal bovine serum, 
100 U/ml of penicillin and 100 ug/ml of streptomycin (Life 
Technologies). For the AGS 3D culture, flasks were 
coated with 50ul/cm 2 of matrigel basement membrane 
(BD Biosciences) at a 1.5:1 proportion to serum-free 
medium. A 2 x 10 4 cells/cm 2 suspension was seeded on 
top and maintained for 14 days with medium change 
every 2 days. For the Caco-2 3D culture, coverslips were 
coated with 60ul/cm 2 of matrigel. A 6 x 10 3 cells/cm 2 
suspension plus 2% matrigel was seeded on top and main- 
tained for 8 days with medium change every 2 days. 
To quantify lumen formation, > 100 cysts were microscop- 
ically examined. For proteasome inhibition, cells were 
treated with 25 uM MG132 (Calbiochem) or vehicle 
treated with DMSO. To inhibit transcription, cells were 
exposed to 10Lig/ml of Actinomycin D (Sigma). 

Constructs and site-directed mutagenesis 

The previously published pCMV-MEX3A expression 
vector (28) was used together with a pCMV-Tag3B 
empty vector (Agilent Technologies) in transfections. 
A pRLControl construct containing a humanized Renilla 
luciferase (Rluc) coding sequence (35) was used as a 
backbone to create the pRLCDX2 vector, encoding a 
luciferase fusion transcript to the parental CDX2 
3'untranslated region (UTR). The QuickChange site- 
directed mutagenesis kit (Stratagene) was then used to 
introduce specific mutations in the previous plasmid, to 
generate the pRLACDX2 construct, with a mutated 
MEX-3 recognition element (MRE). Oligonucleotides 
containing the desired mutations were designed according 
to the manufacturer's instructions (Supplementary Table 
SI). Taking into account the degenerate consensus 
sequence described for MEX-3 binding in C. elegans, a 
mutational background of Cytidine was used because pre- 
sumably this base is not tolerated in the MRE (36). 
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Transfections, RNAi and luciferase assays 

Transient transfections were done using Lipofectamine 
2000 reagent according to the manufacturer's guidelines 
(Life Technologies). A DNA (ug) to Lipofectamine 2000 
reagent (ul) ratio of 1:1.5 in OPTI-MEM medium was 
used for routine transfection experiments of expression 
vectors. For stable transfections, selection was initiated 
48 h post-transfection in medium supplemented with 
0.6mg/ml of G418 (Sigma). Neomycin-resistant positive 
clones obtained through limiting dilution were routinely 
maintained with 0.2mg/ml of G418. A commercial set of 
three Stealth Select siRNA duplexes (HSS 150674, 
HSS 150675 and HSS 150676) directed against human 
MEX3A (Life Technologies) and a custom set of three 
siRNA duplexes directed against CDX2 were used with 
scrambled controls. An siRNA duplexes (pmol) to 
Lipofectamine 2000 reagent (ul) ratio of 20:1 in OPTI- 
MEM medium was used for inhibition experiments. 
Luciferase reporter assays were performed with the Rluc 
Assay System (Promega), and (3-galactosidase activity was 
used for normalization of experimental variations. 

RNA isolation and quantitative real-time PCR 

Total RNA was extracted using TRI Reagent (Sigma) and 
reverse transcribed using the Superscript II Reverse 
Transcriptase kit (Life Technologies). Analysis of BMI1, 
CDX2, GAPDH, LGR5, MEX3A and Rluc mRNA 
expression was performed in an ABI Prism 7500 system 
using SYBRgreen reagent (Life Technologies) and specific 
primer pairs (Supplementary Table SI). Each sample was 
amplified in triplicate and specificity confirmed by dissoci- 
ation analysis. The 18S rRNA expression was measured 
for normalization of target gene abundance. 

Microarrays and data processing 

Experiments were performed at the National Facility for 
DNA Microarrays (Aveiro, Portugal). Three independent 
2D and 3D AGS cell cultures were selected for total RNA 
extraction using TRI Reagent. RNA quantity and quality 
were assessed using the Nanodrop ND-1000 (Thermo 
Scientific) and 2100 Bioanalyzer (Agilent Technologies) 
systems, and only samples with a RNA integrity number 
above nine were considered for further study. Preparation 
and labelling of the RNA was performed using the 
One-Color Microarray-Based Gene Expression Analysis 
Quick Amp Labelling kit (Agilent Technologies). Briefly, 
600 ng of total RNA from 2D and 3D biological replicas 
was used as input together with spike-in controls to 
generate Cyanine 3-labelled cRNA. The amplified cRNA 
samples were purified with RNeasy Mini kit (Qiagen) and 
hybridized to Human Gene expression 4 x 44K v2 
Microarray slides (Agilent Technologies) at 65°C for 
17 h. After washing procedures, images of the hybridiza- 
tions were acquired using a G2505B Microarray Scanner 
(Agilent Technologies). The Feature Extraction software 
was used for spot identification, background subtraction 
and quantification of the fluorescent signal. Raw expres- 
sion data were processed using BRB-ArrayTools v3.8.1 
software (37). After base 2 log transformation and 



average of probe replicates, normalization was applied 
using the median intensity over the entire array to 
minimize systematic variance. Differentially expressed 
transcripts between the two culture conditions were 
filtered using the Class Comparison tool, performing an 
unpaired sample /-test with a P<0.01 and considering a 
minimal 1.5-fold change. Hierarchical clustering analysis 
of significantly altered genes was conducted in TIGR 
MultiExperiment Viewer v4.8.1 software (38) using 
Euclidean correlation and average linkage clustering, and 
expression values indicated colorimetrically. Functional 
annotation of differentially expressed genes, identified by 
extending the unpaired /-test parameter to P < 0.05 and 
maintaining the minimal fold change was performed 
with DAVID program (39). Microarray data have been 
submitted to the ArrayExpress database (http://www.ebi. 
ac.uk/arrayexpress/) and assigned the identifier E-MTAB- 
1234. The complete list of differentially expressed genes is 
detailed in Supplementary Table S2. 

Protein extraction and western blot analysis 

Cells were lysed for 30 min on ice in a lysis buffer 
containing 20 mM Tris-HCl (pH 7.5), 150mM NaCl, 
2mM ethylenediaminetetraacetic acid and 1% IGEPAL 
(Sigma), supplemented with Complete protease inhibitor 
cocktail (Roche Applied Science), 1 mM PMSF and 1 mM 
Na 3 V0 4 . Lysates were centrifuged at 12 000 rpm for 20 
min at 4°C and the supernatant recovered. Protein con- 
centration was determined using the BCA Protein Assay 
Reagent (Thermo Scientific). Protein extracts (30-50 ug) 
were run on 10% sodium dodecyl sulphate-polyacryl- 
amide gel electrophoresis, transferred to a nitrocellulose 
membrane, blotted overnight with appropriate antibodies 
(Supplementary Table SI), and signals revealed with 
ECL detection kit (GE Healthcare Life Sciences). Actin 
levels were used to normalize protein expression, and 
quantification of western blots was performed using Fiji 
software (40). 

Flow cytometry 

Caco-2 cells were harvested 48 h after transfection at con- 
fluence and fixed with 1% paraformaldehyde for 20 min 
at room temperature followed by permeabilization with 
0.1% Triton X-100 (Sigma) for 5 min on ice. Staining 
was performed with anti-MEX3A antibody and visualized 
with goat anti-rabbit FITC-conjugated secondary 
antibody. For DNA content assessment, cells were 
incubated with a propidium iodide (50 ug/ml) and 
RNAse A (200 ug/ml) solution (Sigma) for 30 min at 
room temperature. Samples were read in a FACSCanto 
II (BD Biosciences) flow cytometer, and analysis was 
performed using the FlowJo software. 

RNA-immunoprecipitation assay 

Cells were lysed 48 h after transfection as indicated for 
protein extraction but with the addition of 20 U/ml of 
RNAseOUT Ribonuclease inhibitor (Life Technologies). 
Before the lysis procedure, cells were washed with phos- 
phate buffered saline (PBS) and subjected to ultraviolet 
crosslink (254-nm wavelength) at an energy level 
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of 1.8J/cm 2 on ice for 20 min. Immimoprecipitation 
from protein extracts (750 ug) was performed overnight 
with anti-myc antibody (Clontech) or normal mouse 
immunoglobulin G (IgG) (Millipore) and RNA- 
immunocomplexes recovered by incubation with protein 
G-sepharose bead slurry (Sigma) for 2h in the following 
day. After five washes in lysis buffer containing 
RNAseOUT, RNA was extracted from the bead slurry 
as indicated before and resuspended in 10 ul of dH 2 0. 
RT-quantitative real-time PCR (qPCR) analysis was 
directly performed on the RNA as described. Target 
mRNA levels immunoprecipitated in the MEX3A trans- 
fected cells were normalized to the levels recovered in the 
empty vector transfected cells. 

Immunofluorescence and microscopy 

Cells grown on coverslips were fixed using ice-cold 
methanol for 10 min. For Caco-2 3D assay, cysts were 
fixed with a 2% paraformaldehyde and 0.05% glutar- 
aldehyde solution at room temperature for 20 min, 
followed by quenching with 0.1% NaBH 4 and permeabil- 
ization with 0.5% Triton X-100. Paraffin-embedded 
wild-type murine tissue sections (4 um) were deparaffinized 
with clear-rite (Thermo Scientific) and rehydrated. All 
samples were blocked with non-immune serum diluted in 
10% BSA-PBS, followed by incubation with appropriate 



antibodies (Supplementary Table SI) diluted in 5% BSA- 
PBS. Nuclei were counterstained with 4',6-diamidino- 
2-phenylindole (DAPI, Sigma), and coverslips mounted 
on microscope slides with Vectashield (Vector). Imaging 
acquisition was based on optical sectioning and adjusted 
for brightness and contrast with Fiji software (40). Imaging 
was performed on an Axio Imager Zl (Zeiss) Fluorescence 
microscope with an ApoTome attachment (Figures 3, 7 
and 8 and Supplementary Figure S4), and on a TCS 
SP5 II (Leica) Laser Scanning Confocal microscope 
(Supplementary Figure SI). The middle focal plane of 
Z-stack images with similar total number of planes were 
used for representation of cysts polarity. 

Statistics 

Each experiment was carried out in triplicates at least 
three times, and data were expressed as means ± SD. 
Statistical analysis was performed using Student's f-test. 
A P < 0.05 was considered as significantly different. 

RESULTS 

CDX2 expression is downregulated in AGS cells 
cultured in matrigel 

In an effort to study CDX2 regulation in a physiologically 
relevant microenvironment, we compared CDX2 
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Figure 1. Establishment of AGS 3D cultures and characterization of CDX2 expression. (A) Schematic representation of the culture systems and 
cellular morphology at culture day 4 in bright field microscopy (original magnification, xlOO). (B) Western blot and qPCR of CDX2 expression in 
2D and 3D after 2 weeks culture. Cardinal numbers represent biological replicates. Values for CDX2 mRNA expression in 2D culture were referred 
to as 1. (C) Western blot of CDX2 protein expression on treatment with MG132. (D) qPCR of CDX2 mRNA expression on treatment with 
actinomycin D. Expression levels in the absence of treatment were set at 100%. Depicted half lives were calculated using exponential regression 
(**P = 0.003). 
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expression in AGS cells cultured in a standard 2D condi- 
tion versus a 3D condition with matrigel. Cells depicted 
distinct growth patterns amongst both settings, forming a 
compact monolayer and conspicuous multicellular aggre- 
gates, respectively (Figure 1A). Although there was a 
dramatic decrease in CDX2 protein levels in the 3D 
culture, mRNA showed heterogeneous levels without cor- 
relation with protein (Figure IB). Treating cells with the 
proteasome inhibitor MG132 indicated that CDX2 
downregulation was not related to an increase in 
proteasome-dependent degradation (Figure 1C). On the 
other hand, we found a significant difference in CDX2 
mRNA stability after actinomycin D treatment to 
inhibit transcription (Figure ID). The CDX2 mRNA 
half-life value was higher for the 2D culture 
(166.9 ± 15.5 min) than for the 3D culture (62.9 ± 22.1 
min). These results suggest the existence of a post- 
transcriptional regulatory mechanism, probably acting 
over CDX2 mRNA. 

With the objective of pinpointing the mechanism 
underlying CDX2 downregulation, we performed a 
whole-genome expression array, comparing the transcrip- 
tome of AGS cells in the two culture conditions. Among 
the 340 differentially expressed transcripts, a 1.7-fold 
upregulation of MEX3A was observed in the 3D condition 
(Figure 2A and Supplementary Table S2). Additionally, 
functional annotation analysis showed alterations related 
to fundamental cellular processes like proliferation, differ- 
entiation and survival, namely an increase in cell cycle 



elements, Notch and Wnt pathways and a decrease in 
MAPK and Jak-STAT signalling in the 3D model 
(Figure 2B). An increase in RhoGTPases and integrin 
signalling was also registered, supporting the occurrence 
of increased cell-matrix interactions. Changes in metabolic 
patterns were evident by an increase in amino acid biosyn- 
thesis (Figure 2B). MEX3A is one of the human ortho- 
logues of MEX-3, an RBP with translational inhibitory 
function known to target the CDX orthologue pal-1 
during C. elegans embryonic development. A novel inter- 
action between human MEX3A and CDX2 in a gastro- 
intestinal context could constitute a plausible explanation 
for the reduction of CDX2 protein levels observed in 3D 
culture. Accordingly, MEX3A expression was validated by 
qPCR (Figure 2C). 

MEX3A overexpression leads to CDX2 downregulation 
in AGS cells 

To define whether CDX2 expression is regulated by 
MEX3A, AGS cells cultured in standard conditions were 
transiently transfected with a myc-tagged MEX3A expres- 
sion vector and collected 24 and 48 h later. We showed that 
the myc-tagged MEX3A fusion protein was detectable on 
transfection and accompanied by a 77 and 64% decrease 
in CDX2 protein expression, respectively (Figure 3A). 
In contrast, CDX2 mRNA expression did not change 
(Figure 3B). This result was confirmed in AGS cells 
stably transfected with the myc-tagged MEX3A construct 
up until 96 h of culture (Figure 3C and D). 
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Figure 2. Transcriptome profiling of 2D and 3D AGS cultures. (A) Microarray heat map displaying differentially expressed transcripts (P<0.01, 
fold change > 1.5). Values obtained for MEX3A are highlighted. (B) Functional annotation analysis depicting the preponderance of differentially 
expressed genes in terms of signaling networks (P<0.05, fold change > 1.5). (C) qPCR validation of MEX3A expression. Values for MEX3A mRNA 
expression in 2D culture were referred to as 1. 
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Figure 3. CDX2 regulation by MEX3A in AGS cells. (A) Western blot of transient transfections with a myc-tagged MEX3A expression vector. 
(B) qPCR of CDX2 and MEX3A mRNA expression in the previous samples. The value for CDX2 mRNA expression in the empty vector transfected 
cells at 24 h was referred to as 1. (C) Immunofluorescence for MEX3A and CDX2 in MEX3A stably transfected cells at 48 h culture (original 
magnification, x 400; scale bar 20 urn). (D) Western blot of MEX3A and CDX2 expression in stably transfected cells at different time points. 



It was previously published that MEX3A and MEX3B 
were novel components of P bodies (28,29), discrete cyto- 
solic foci characterized in yeast and mammals as centres for 
translational silencing and mRNA decay (41). In our 
setting, on MEX3A-induced CDX2 protein decrease, we 
observed a specific subcellular localization of MEX3A in P 
bodies, as confocal microscopy analysis revealed a partial 
overlap between myc-tag staining and two known P body 
components, the human mRNA-decapping enzyme 1A 
(DCP1A) and the human enhancer of mRNA-decapping 
protein 4 (EDC4) (Supplementary Figure SI). Together, 
these results demonstrate that MEX3A has the ability to 
regulate CDX2 expression at the post-transcriptional level, 
probably by interplay with P bodies. 

MEX3A interacts with a MRE present in CDX2 3'UTR 

A bioinformatics search conducted throughout the CDX2 
3'UTR led to the identification of the MEX-3 degenerate 
consensus binding sequence (36). The MRE present in 



CDX2 mRNA was defined as a bipartite element that 
consists of AGAG and UUUA motifs separated by two 
Uracil bases (Figure 4A). To assess whether MEX3A 
associates with CDX2 3'UTR, we preserved RNA 
protein interactions by ultraviolet covalent linkage in 
AGS cells transiently transfected with the MEX3A expres- 
sion vector and performed an RNA-immunoprecipitation 
(RIP) assay. We found a significant enrichment in CDX2 
mRNA recovered with the anti-myc antibody compared 
with control IgG (Figure 4B). On the contrary, no enrich- 
ment for GAPDH mRNA was observed (Figure 4B). 
Likewise, no enrichment for TBP mRNA and 18S 
rRNA was detected (data not shown). Next, we asked 
whether MEX3A binds specifically to the MRE present 
within the CDX2 3'UTR. To address this question, Rluc 
expression vectors containing the wild-type CDX2 3'UTR 
(pRLCDX2) or the CDX2 3'UTR with a mutated MRE 
(pRLACDX2) were used in an RIP assay as before 
(Figure 4C). A complete loss in MEX3A binding 
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Figure 4. Mechanistics of MEX3A interaction with CDX2 mRNA in AGS cells. (A) Terminal sequence of CDX2 mRNA, with the predicted MRE 
in the 3'UTR highlighted in bold and underlined. (B) qPCR showing CDX2 and GAPDH mRNAs immunoprecipitated with myc-tag antibody in 
MEX3A transfected cells (*P = 0.018). The values for CDX2 and GAPDH mRNA levels in the IgG sample were referred to as 1. (C) Schematic 
representation of pRL constructs. (D) qPCR showing Rluc mRNA immunoprecipitated with myc-tag antibody for pRLCDX2 transfection 
(****/>< 0.0001). Values for Rluc mRNA expression in IgG samples were referred to as 1. (E) Luciferase activity assay for the different pRL 
constructs (**P = 0.006 for pRLControl/pRLCDX2 and *P = 0.036 for pRLCDX2/pRLACDX2). The values obtained for luciferase expression in 
the pRLControl co-transfected cells were referred to as 1. 



capacity was detected when the MRE was mutated, as 
observed by comparing the levels of Rluc mRNA re- 
covered with anti-myc antibody between the pRLCDX2 
and pRLACDX2 transfected cells (Figure 4D). 
Additionally, luciferase activity assays proved that the 
construct with the mutated MRE was insensitive to 
MEX3A function, as a significant reduction in Renilla ex- 
pression was only detected with the pRLCDX2 transfec- 
tion (Figure 4E). We conclude that MEX3A is able to 
interact with CDX2 mRNA and that the major determin- 
ant for this interaction is a putative MRE located in the 
transcript 3'UTR. 

Intestinal phenotype of Caco-2 cells is affected by 
modulating MEX3A levels 

To assess the functional consequences of the previous ob- 
servations and confirm that they were cell-type independ- 
ent, we studied CDX2 regulation by MEX3A in another 
model. Caco-2 cell line was chosen for its ability to spon- 
taneously differentiate into an enterocytic-like phenotype 
on reaching confluence, which has turned it into the most 
widely used system to study intestinal cell maturation (42). 



Furthermore, this process was previously shown to be 
accompanied by an increase in CDX2 protein expression 
that, in some established Caco-2 clones, is not explained 
by higher transcriptional levels (20). First, we evaluated 
the timing of CDX2 and MEX3A expression in Caco-2, 
harvesting cells at pre-confluence (day —2), confluence 
(day 0) and post-confluence (day 2-8) time points. 
Consistent with previous reports, CDX2 protein expres- 
sion progressively increased starting from confluence and 
during differentiation (Supplementary Figure S2A). This 
trend was also verified at the mRNA level, though not to 
the same extent (Supplementary Figure S2B). By contrast, 
MEX3A mRNA expression varied in a complementary 
manner to that of CDX2, with a decreased expression 
associated to more differentiated cells (Supplementary 
Figure S2B), although a protein decrease was only 
visible at day 8 (Supplementary Figure S2A). We inhibited 
endogenous MEX3A using siRNA and observed that 
CDX2 protein started to be highly expressed at earlier 
time points in cells with downregulated MEX3A expres- 
sion (Figure 5A). Paradoxically, this effect was reverted 
from confluence onwards, with lower CDX2 protein levels 
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detected at later time points in cells with MEX3A inhib- 
ition (Figure 5A). To evaluate whether this result was 
dependent on MEX3A activity or due to transcriptional 
regulation, CDX2 mRNA expression was assessed. As 
expected, CDX2 mRNA levels were not altered in the 
pre-confluence time points but were downregulated from 
confluence onwards (Figure 5B), suggesting that after con- 
fluence, other MEX3A regulatory functions indirectly 
interfere with endogenous CDX2 levels. 

It is known that Caco-2 cells acquire expression of dif- 
ferent intestinal columnar lineage markers on differenti- 
ation, such as Villin. To determine the effect of MEX3A 
overexpression in intestinal differentiation, we generated a 
Caco-2 cell line stably transfected with the myc-tagged 
MEX3A construct. Expression of the myc-tag was con- 
firmed in different confluence states, and a concomitant 
reduction in CDX2 protein levels was observed 
(Figure 5C). Furthermore, we detected decreased Villin 
expression in MEX3A overexpressing cells (Figure 5C). 



As differentiation and growth arrest are spatially and tem- 
porally related events in the intestinal epithelium, we 
studied the expression of cyclin Dl, a positive regulator 
of proliferation required for cell cycle Gl/S transition. 
Sustained levels of cyclin Dl were detected in MEX3A 
overexpressing cells in latter time points, when the 
culture is normally already less proliferative, as observed 
by the striking decrease of cyclin Dl in mock cells 
(Figure 5C). To ascertain the cell cycle distribution 
related to MEX3A overexpression, we transiently trans- 
fected near-confluent Caco-2 cells with the myc-tagged 
MEX3A construct and performed DNA content analysis 
by flow cytometry. Measurements indicated that cells 
overexpressing MEX3A possess a different cell cycle 
profile (Figure 5D), with a reduced G0/G1 population 
(32 against 49% in the empty transfected cells) and an 
increased S phase population (41 against 30% in the 
empty transfected cells). Given the observed alterations 
in differentiation and proliferation parameters, we 
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Figure 5. MEX3A modulation of intestinal phenotype in Caco-2 cells. (A) Western blot of MEX3A inhibition with transfection of siRNAs per- 
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cycle phases of the empty vector transfected cells and MEX3A-transfected cells (gate MEX3A exogenous) are shown below (*P = 0.03 for GO/ 
Gl phase and *P = 0.05 for S phase). 
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sought to determine whether MEX3A overexpressing cells 
possessed intestinal 'progenitor or stem-like' features. 
Assessment of expression of well-established intestinal 
stem cell markers (43^16) revealed that these cells show 
overall increased levels of LGR5, BMII and MSI1 
(Figure 6), whereas OLFM4 basal expression was barely 
detectable and did not change (data not shown). 
Importantly, this seems to be independent of CDX2, as 
siRNA-mediated CDX2 downregulation did not repro- 
duce the increase in MSI1 levels (Figure 6C). 

Caco-2 cells retain the ability to polarize and form a 
transporting epithelial monolayer in culture (42). This po- 
larization process is intimately connected with the forma- 
tion of cell-cell contacts and adhesion. In standard culture 
conditions, Caco-2 cells overexpressing MEX3A showed 
alterations in the distribution of the tight junction marker 
ZO-1 compared with the mock cells (Figure 7A). We then 
used a 3D cell cyst assay (4) in matrigel to compare the 
phenotypic response of both cell lines. We observed that 
mock cells developed cyst-like structures in a few days, 
with ~50% showing a well-defined central lumen 
(Figure 7B-D). In contrast, >80% of the MEX3A- 
expressing cysts, in which CDX2 reduction and MEX3A 
increase was confirmed, failed to elaborate such a lumen 
(Figure 7B-D). Additionally, we observed distinct 
patterns of E-cadherin and Phalloidin expression. We 
detected basolateral E-cadherin expression and accumula- 
tion of filamentous-actin (F-actin) around the hollow 
lumen in the mock cysts (Figure 7D), and a reduced 
E-cadherin expression in the MEX3A-expressing cysts, 



with F-actin appearing ubiquitously distributed (Figure 
7D). Thus, we show that the establishment of proper 
apical-basal identity is affected on MEX3A 
overexpression. 

MEX3A is differentially expressed in the crypt-villus unit 

To address MEX3A relevance in vivo, we analysed its 
presence in mouse normal intestine and compared it 
with CDX2 expression. We observed that MEX3A is 
expressed in small intestine and colon, exhibiting mainly 
a nuclear staining pattern. It shows a stronger expression 
in the lower portion of the crypt-villus unit in small intes- 
tine and reduced expression in the uppermost part of the 
colonic crypts (Figure 8). In the small intestine, CDX2 is 
expressed in the same cells as MEX3A except in the base 
of the crypts (Figure 8, insert) that exhibit mainly MEX3A 
protein. In colon, CDX2 is more expressed in the surface 
cells (Figure 8). These results suggest that MEX3A expres- 
sion might be important to define intestinal differentiation 
patterns. 

DISCUSSION 

In this study, we provide the first demonstration that 
MEX3A represses CDX2 expression in the gastrointes- 
tinal context, putatively as a translational repressor, with 
direct implications in intestinal differentiation, polarity 
and sternness features, key mechanisms for tissue homeo- 
stasis that are frequently altered in pathological conditions 
and inextricably linked with cancer. 



Nucleic Acids Research, 2013, Vol. 41, No. 7 3995 




MEX3A 



f m 



B 



48h 



72h 9 days 



■ Single central lumen I No single central lumen 

100 -i 









CDX2 / MEX3A / DAP 


3D cysts 

CO 00 

o o 








&■ -3. - ' i 


«J 

2 40 - 
o 

20 ■ 




Mock MEX3A 
Caco-2 cell line 



Mock 



MEX3A 



Phalloidin 


Merge 






E-cadherin 

■■ 

1 




DAPI 









Phalloidin 


Merge 


. " i -*X _ (* I 








E-cadherin 






; - V -<? "•-*•/•' 


DAPI 


\ - 4-* V 







Figure 7. Polarity alterations induced by MEX3A in Caco-2 cells. (A) Immunofluorescence showing ZO-1 expression in Caco-2 mock and MEX3A 
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3D cell culture systems have been shown to enable 
physiological and functional differentiation of several epi- 
thelial cell types (47,48), constituting a promising alterna- 
tive to overcome standard cell culture limitations. 
Accordingly, our transcriptomic analysis performed on a 
3D model of AGS cells, which have a significant level of 
endogenous CDX2 responsive to different molecular 
stimuli (13-15), allowed us to disclose MEX3A as a mo- 
lecular player involved in the regulation of CDX2 transla- 
tion. The array presented an increased MEX3A expression 
in 3D culture and no alteration in the levels of the other 



MEX3 family members. We further demonstrated that 
MEX3A overexpression leads to marked CDX2 protein 
decrease in two cell lines. We tried to overcome the limi- 
tation of using overexpression systems by using two dif- 
ferent cell lines that were both transiently and stably 
transfected with a MEX3A plasmid, which gave concord- 
ant results. In addition, CDX2 negative regulation by 
MEX3A was confirmed in a more physiological context 
using a siRNA approach towards endogenous MEX3A in 
Caco-2 cells. We proved that MEX3A is able to interact 
with CDX2 mRNA through a canonical MRE present 
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in the 3'UTR. This sequence, which we now show to be 
functionally relevant in humans, seems to be the single 
determinant of MEX3A binding, independently of the 
upstream 5'coding region. In fact, bioinformatics 
analysis shows that the degenerate MRE motif has been 
evolutionarily conserved in different CDX2 homologues 
(Supplementary Figure S3), namely in chimpanzee, 



mouse, rat, zebrafish, fruit fly and frog, suggesting that 
MEX3A is critical for CDX2 regulation. 

Endogenously, MEX3A was mainly localized in the 
nuclear compartment of intestinal epithelial cells, both 
in vitro and in vivo. This result is crucial to show that 
a biological background for the regulation of CDX2 by 
MEX3A exists. In this regard, a variant form of MEX3D 
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called TINO, which has been shown to negatively regulate 
BCL-2 expression by transcript destabilization, is also pre- 
dominantly localized in the nuclei of HeLa cells (33). 
On the other hand, MEX3A transfectants showed pre- 
dominant cytoplasmic staining. A potential problem 
could be antibody specificity; however, this was success- 
fully evaluated by transfection of siRNA duplexes directed 
against MEX3A (Supplementary Figure S4). Therefore, 
the differential localization in distinct expression 
backgrounds might be because a certain threshold in ex- 
pression levels has to be achieved for cytoplasmic trans- 
location to occur, given that MEX3 proteins are capable 
of performing nucleocytoplasmic shuttling (28). On the 
other hand, specific subcellular distribution might be 
related with protein phosphorylation, in accordance with 
published data describing MEX3 members as phospho- 
proteins (28,29). Whichever the case, this does not contra- 
dict post-transcriptional regulation, as it can be elicited at 
multiple points of the transcript lifespan, including 
pre-mRNA processing in the nucleus, export from the 
nucleus to the cytoplasm and subsequent coordinated 
trafficking of the mature mRNA to the translation ma- 
chinery. It also remains to be fully clarified the conse- 
quence of partial accumulation in P bodies, structures 
involved in processes of mRNA degradation, nonsense- 
mediated mRNA decay, translational repression and 
RNA-mediated gene silencing (49), although there are 
several RBPs with established roles in translational regu- 
lation known to colocalize with P bodies (50,51). 

We used the Caco-2 cell line model to modulate MEX3A 
expression and assess its phenotype. Surprisingly, MEX3A 
inhibition produced distinct effects over CDX2 levels de- 
pending on cellular confluence. MEX3A might selectively 
regulate unique subsets of targets in different culture con- 
ditions. In agreement, it is known that transcriptomic and 
proteomic changes occur during the progression from the 
proliferative state to spontaneous differentiation of Caco-2 
cells. Conversely, MEX3A might have a dual role, acting 
both as a repressor or enhancer contingent on the cellular 
microenvironment, as observed for other RBPs, like HuR 
(52). It is possible that intricate associations of MEX3A 
with other molecular effectors, namely other proteins or 
microRNAs, determine divergent regulatory endings. 
MEX3A overexpression in Caco-2 cells resulted in 
pronounced phenotypic alterations. CDX2 and Villin 
downregulation were indicative of a loss in intestinal dif- 
ferentiation, which was further confirmed by the flow 
cytometric profile, showing reduced G0/G1 population, 
usually associated with less-differentiated cells. Another 
hallmark feature of MEX3A overexpression was the 
altered cellular polarization in standard culture, as well 
as the impaired ability to form polarized structures in the 
presence of matrigel. These effects may be mediated by 
CDX2, as this transcription factor was previously shown 
to regulate intestinal Villin through recruitment of the 
Brm-type SWI/SNF complex to its promoter (53), and 
MEX3A-expressing cysts closely resemble the ones 
obtained with Caco-2 cells in which CDX2 suppression 
was achieved by lentiviral short-hairpin RNA particles 
(4). It is, therefore, important to ascertain the biological 
setting where this regulation might be determinant, which 



is suggested by the predominant expression of MEX3A in 
the stem, transit-amplifying and migrating post-mitotic 
cells of the intestine. Although CDX2 protein can be 
detected in most of these cells, its level is lower in the 
crypts compared with the uppermost differentiated cells 
of the villi (2,19,20). A similar increasing bottom-up 
gradient has been described along the colonic gland axis 
in the distal colon epithelium. By contrast, in situ hybrid- 
ization revealed that CDX2 mRNA was homogeneously 
distributed along the entire crypt-villus axis (20). Given 
this CDX2 protein expression gradient, that we also 
show, and lack of correlation with mRNA, it is likely 
that MEX3A fine-tunes CDX2 levels in vivo as well, in a 
transcription-independent manner, providing swift avail- 
ability of the protein to meet the physiological require- 
ments of the continuously renewed gut epithelium. 
Moreover, it has been shown that intestinal stem cells 
cannot differentiate into any of the intestinal lineages in a 
background of Cdx2 ablation (54), revealing the need for 
tight CDX2 regulation as determinant to proper phenotype 
switching. Most interestingly, we observed higher expres- 
sion of different intestinal stem cell markers when we 
overexpressed MEX3A in Caco-2, suggesting that this 
protein is associated with stem cell features. Our AGS 3D 
model also showed a significant upregulation of OLFM4 
(Supplementary Table S2), previously identified as a 
marker for LGR5+ stem cells in human intestine (46), 
although this gene was not directly upregulated in Caco-2 
cells, which might be due to intrinsic properties of each cell 
line. Strengthening the hypothesis for a role of MEX3A in 
stem cell potential, a recent publication showed that 
MEX3A is part of the molecular signature of the LGR5+ 
intestinal stem cells, presenting a 1.64-fold increase in 
relation to daughter cells, along with MSI1, LGR5 and 
OLFM4 (55). Furthermore, our study reinforces the 
increasing knowledge that other regulatory mechanisms 
in addition to transcriptional ones have important func- 
tions in stem cells. Although we do not know yet how 
MEX3A relates with the intestinal stem cell phenotype, 
we hypothesize that MEX3A overexpression in itself, 
together with the induction of a lessened differentiated 
phenotype and polarity defects mediated by CDX2 
downregulation, might be critical to allow a permissive en- 
vironment for the appearance of sternness features. 
Furthermore, the MEX3A-CDX2 axis might be important 
during early embryogenesis where CDX2 is required for 
correct trophectoderm differentiation, while absent from 
the inner cell mass, in a process that requires tight regula- 
tion (56). Interestingly, the pattern of expression of the C. 
elegans orthologues of both proteins was found to be 
mutually exclusive in early stages of embryogenesis 
(24,27). Still, it remains to be assessed the relevance of 
MEX3A in multiple pathological contexts of the gastro- 
intestinal tract where differentiation abnormalities 
directed by CDX2 are key events. 

In conclusion, we have identified a novel role for 
MEX3A protein in the regulation of intestinal differenti- 
ation, polarity and sternness features, partially mediated 
by the repression of CDX2. This is the first description of 
a CDX2 regulatory mechanism based on its mRNA 
control by an RBP, having a significant impact in 
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intestinal homeostasis and likely in gastrointestinal 
carcinogenesis. 
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